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ABSTRACT 

We measure the matter power spectrum from 31 Lya spectra spanning the redshift 
range of 1.6-3.6. The optical depth, r, for Lya absorption of the intergalactic medium 
is obtained from the flux using the inversion method of Nusser & Haehnelt (1999). 
The optical depth is converted to density by using a simple power law relation, r oc 
(1 + S) a . The non- linear ID power spectrum of the gas density is then inferred with 
a method that makes simultaneous use of the 1 and 2 point statistics of the flux and 
compared against theoretical models with a likelihood analysis. A Cold Dark Matter 
(CDM) model with standard cosmological parameters fits the data well. The power 
spectrum amplitude is measured to be (assuming a flat Universe), as = (0.9 ± 0.09) x 
(f2 m /0.3) -0 ' 3 , with a varying in the range of 1.56 — 1.8 with redshift. Enforcing the 
same cosmological parameters in all four redshift bins, the likelihood analysis suggests 
some evolution in the density-temperature relation and the thermal smoothing length 
of the gas. The inferred evolution is consistent with that expected if reionization of 
Hell occurred at z ~ 3.2. A joint analysis with the WMAP results together with a 
prior on the Hubble constant as suggested by the HST key project data, yields values 
of f2 m and as that are consistent with the cosmological concordance model. We also 
perform a further inversion to obtain the linear 3D power spectrum of the matter 
density fluctuations. 

Key words: cosmology: theory - intergalactic medium - hydrodynamics - large-scale 
structure of Universe - quasars: absorption lines 



1 INTRODUCTION 

The numerous Lya absorption features observed in quasar 
spectra blue-wards of their Lya emission line known as the 
Lya forest, provide one of the main probes of the intergalac- 
tic medium (hereafter IGM) (Bahcall & Salpeter 1965; Gunn 
& Peterson 1965). In recent years two main advances have 
shaped the accepted view on the origin of the Lya absorbing 
structures. First, the advent of 10-meter class telescopes 
equipped with high-resolution echelle spectrographs (HIRES 
on Keck and UVES on the Very Large Telescope) has pro- 
vided us with data of unprecedented quality (see Rauch 
1998 for a review). Second, the emergence of a theoretical 
paradigm within the context of the cold dark matter (CDM) 
cosmology (e.g. Bi, Boerner & Chu 1992) supported by nu- 
merical hydrodynamical simulations (Cen et al. 1994; Zhang, 
Anninos & Norman 1995; Miralda-Escude et al. 1996; Hern- 
quist et al. 1996; Wadsley & Bond 1996; Zhang et al. 1997; 
Theuns et al. 1998; Machacek et al. 2000; Viel, Haehnelt & 
Springel 2004; Tytler et al. 2004; Jena et al. 2005) and semi- 



analytical studies (e.g. Pichon et al. 2001; Matarrese & Mo- 
hayaee 2002; Viel et al. 2002). According to this paradigm, 
the absorption is produced by volume filling photoionized 
gas that contains most of the baryons at redshifts z ~ 3 (see 
e.g. Efstathiou, Schaye & Theuns (2000) for a recent review) , 
where the absorbers are locally overdense extended struc- 
tures, close to local hydrostatic equilibrium (Schaye 2001). 
The paradigm also predicts that most of the gas probed 
by the Lya forest - absorption features with column den- 
sity ^ 10 13 ' 5 cm -2 - resides in mildly non-linear dark-matter 
overdensities. 

On scales smaller than the Jeans scale the baryonic gas 
is smoothed by pressure forces, erasing the small-scale fluc- 
tuations of the gas density, and setting its distribution apart 
from the dark matter component. On these small scales 
the width of the absorption features is determined by the 
gas thermodynamical properties enabling the measurement 
of the IGM temperature and temperature-density relation 
(e.g., Schaye et al. 2000, Theuns & Zaroubi 2000, McDon- 
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aid et al. 2001, Theuns et al. 2002a and 2002b, Gleser et 
al. 2005). 

On scales larger than the Jeans scale, however, the gas 
distribution faithfully follows that of the underlying dark 
matter. The gas distribution on these scales provides a probe 
of the dark matter distribution and its power spectrum 
(Croft et al. 1998; Hui 1999; Croft et al. 2002, McDonald et 
al. 2000, Hui et al. 2001; Croft et al. 2002; Viel et al. 2003, 
2004, 2005; McDonald et al. 2005; Seljak et al. 2005; Des- 
jacques & Nusser 2005; Lidz et al. 2005, Viel & Haehnelt 
2005). 

The "standard" method used for measuring the matter 
power spectrum from the Lya forest is based on the result 
obtained in numerical simulations showing that the normal- 
ized flux power spectrum is proportional to that of the un- 
derlying matter (Croft et al. 1998). The calibration of the 
relation between the two relies on numerical simulations and 
on the value of the mean flux. Both of these are somewhat 
uncertain. The calibration of the relation depends on the 
specific cosmological parameters of the simulation, and the 
mean flux on the specific quasar (QSO) spectrum at hand. 

To date, almost all matter power spectra inferred from 
the Lya forest are based on this relation. There is, however, 
some debate on what to adopt for the mean level of the flux. 
Early power spectrum measurements from the Lya forest 
(e.g., Croft et al. 1998; 2002; McDonald et al. 2000) adopted 
a rather low level of the mean flux, and inferred a relatively 
low amplitude of the fluctuations. This result is not com- 
pletely consistent with the large scale angular power spec- 
trum amplitude and the early reionization of the Universe 
inferred from the temperature and polarisation data of the 
WMAP satellite (Kogut et al. 2003). This tension has been 
the primary reason for the WMAP team to suggest a run- 
ning spectral index model (Spergel et al. 2003). A num- 
ber of authors have recently pointed to the strong depen- 
dence of the inferred amplitude on the adopted mean flux 
level and have argued that the errors of dark matter power 
spectrum inferred from Lya forest data have been underesti- 
mated (Zaldarriaga, Scoccimarro & Hui, 2003; Zaldarriaga, 
Hui & Tegmark, 2001; Gnedin & Hamilton 2002; Seljak, Mc- 
Donald & Makarov 2003). 

Recent studies (Kim et al. 2004, Viel et al. 2004 and 
Seljak et al. 2005, Lidz et al. 2005, Viel & Haehnelt 2005) 
have all adopted the higher values of the mean flux sug- 
gested by high-quality absorption spectra. These studies ob- 
tained matter power spectra from Lya data that are consis- 
tent with the WMAP results without the need for the run- 
ning spectral index power spectrum proposed by the WMAP 
team. 

The current study follows a different route in which the 
matter power spectrum is measured by inverting the normal- 
ized Lya flux to obtain the optical depth for Lya absorption 
(Nusser & Haehnelt 1999, 2000; hereafter NH99 and NH00 
respectively). In a system in photoionization- recombination 
equilibrium, such as the IGM, the optical depth for 
Lya absorption is to a good approximation a power law 
function of the underlying density with the power law in- 
dex determined by the temperature-density relation of the 
IGM gas. We use this simple relation to infer the line-of- 
sight distribution of the gas density. The shape of the power 
spectrum and the probability distribution of the gas density, 
are then used to infer the shape and amplitude of the ID- 



power spectrum of the gas density separately, without the 
need for assuming a mean flux level. State-of-the-art hydro- 
dynamical simulations are then used to calibrate and test 
the method. The method is applied to 31 high resolution 
Lya spectra of which the newly acquired LUQAS sample 
constitute the main part. 

The paper is organised as follow: §2 describes the data 
set. §3 shows how we infer the non-linear ID power spec- 
trum of the gas density from the Lya forest data using the 
NH99 & NH00 inversion method. §4 presents the analytical 
modelling of the non-linear ID power spectrum of the gas 
from the linear 3D power spectrum of matter and describes 
the use of numerical hydro-simulations to correct for the ef- 
fect of redshift distortions and gas pressure. The likelihood 
analysis employed to constrain the parameters of the mat- 
ter power spectrum and the thermal state of the IGM, is 
discussed in §5. In §6, our best estimate of the 3D matter 
power spectrum is presented and compared to results from 
previous studies. The main conclusions are given in §7. 



2 THE DATA SET 

The sample used in this study consists of 31 high-resolution 
high S/N spectra. 27 of the spectra were obtained with 
the Ultra- Violet Echelle Spectrograph (UVES) on VLT, 
Paranal, Chile, over the period 1999-2002. The 27 spectra 
were taken from the ESO archive and are publicly avail- 
able (P.I.: J. Bergeron [Bergeron et al. 2004]); this sample 
is known as the LUQAS sample (Kim et al. 2004). The 
LUQAS sample was selected based on the following crite- 
ria: 1) S/N larger than 25 in the Lyman-a forest region; 2) 
complete or nearly complete coverage of the Lyman-a forest 
region; 3) no damped-Lya (DLA) systems in the Lya forest 
region, though few spectra have sub-DLAs (column density 
10 19 ' 0-20 ' 3 cm -2 ); 4) no broad absorption line systems; 5) 
publicly available as of January 1, 2003. The total redshift 
path of the sample is z = 13.75. For more details on the 
LUQAS sample see Kim et al. (2004). The rest of the spectra 
were obtained from various other publicly available spectra 
that fulfill similar criteria. For more details on the remaining 
spectra, see Hu et al. (1995) and Theuns et al. (2002a) and 
references therein. All the 31 spectra used here, have signal- 
to-noise ratios of 40-50 per pixel, and a similar resolution, 
(A/AA > 40000). 

Figure shows the redshift range covered by all the 
spectra used in this analysis. The median redshift of the 
sample is (z) — 2.55 with a cumulative redshift path of about 
16.8. 



3 RECOVERING THE ID POWER 

SPECTRUM OF THE GAS DENSITY IN 
REDSHIFT SPACE 

3.1 From flux to gas density 

Nusser & Haehnelt (NH99 & NH00) introduced an algorithm 
to invert the observed flux in the Lya forest region of QSO 
absorption spectra to obtain the gas density along the line- 
of-sight (see also Pichon et al. 2001). We will use here a 
modified version of this algorithm to measure the ID power 
spectrum of the gas density. This is the first step in our 



© 0000 RAS, MNRAS 000, 000-000 



The Matter Power Spectrum From The Lya Forest 3 



E 

D 
Z 

E 

k. 
o 

<D 
CL 
C/1 



HE1 1 22-1 648 



00453-423 



Q1 101-264 



00329-385 



Q0420-388 



Q0103 



-260 

— ooi 



Q0 1 22-3 80 
Q0302-003 

09-3518 



PKS2126-158 



Q0002-422 
PKS1448-232 



PKS0237-23 



PKS0329-255 



HE2347-4342 



J2233-606 
HE1 347-2457 



HE2217-2818 



HE1 341 



HE0940-1050 



HE1158-1843 



HE0151-4526 

HE0515-4414 



946-7658 



HE00 01 -2340 
HS1 
Q1422+231 

Q0956+122 

AP M08279+52 55 

00636+680 



Q0055-269 



Q0014+813 



1.2 



1.7 



3.7 



4.2 



2.2 2.7 3.2 
redshift range 

Figure 1. The redshift range of the QSO spectra used in this 
paper. 27 of the spectra are taken from the LUQAS QSO sample 
(Bergeron et al. 2004; Kim et al. 2004). QSO 0956+122 is taken 
from Hu et al. (1995). The remaining 4, at the high redshift end 
of the sample, are various other publicly available spectra (see 
Theuns et al. (2002a) and references therein). 



endeavour to constrain the matter power spectrum and the 
thermal state of the IGM. For the sake of completeness, the 
description of the algorithm is repeated here in some detail. 

The optical depth in redshift space due to resonant 
Lya scattering is related to the Hi density along the line 
of sight (LOS) in real space by 



T(Z) = <7 



H{z) 



nm(z, x)TL[z - 



v p (x),b(x)]dx, (1) 



where <to is the effective cross section for resonant line scat- 
tering, H(z) is the Hubble constant at redshift z, x is the 
real space coordinate (in kms -1 ), H is the Voigt profile nor- 
malized such that f H dx = 1, v p (x) is the LOS peculiar 
velocity, and b(x) is the Doppler parameter due to ther- 
mal/turbulent broadening. The absorption features in the 
Lya forest are mainly produced by regions of low to mod- 
erate densities, where photoheating is the dominant heating 
source and shock heating is not important. 

Hydrogen in the IGM is highly ionized (Gunn & Pe- 
terson 1965, Scheuer 1965) and the photoionization equi- 
librium in the expanding IGM establishes a tight correla- 
tion between neutral and total hydrogen density. Numerical 
simulations have supported the existence of this correlation 
and shown that the gas density traces the fluctuations of 
the DM density on scales larger than the Jeans length, so 



that nm = nm 

(«2MM) a . Here nm is the Hi density at 
the mean dark matter density, and the parameter a de- 
pends on the reionization history. The possible range for 
a is 1.56 <C a <C 2 with a value close to 2 just after reion- 
ization, and decreasing at later times (Hui & Gnedin 1997). 
In this relation p DM (x) is the dark matter density smoothed 
on the Jeans length below which thermal pressure becomes 
important. The Jeans length in comoving units in the linear 
regime is given by, 

27TC a 



X 



v^Gp 



1 + z) 



1.3 



0.125 



1.5 



-1/2 



T 



1.5 x 10 4 K 



1/2 



l + (2-a)/0 



1+2 



-1/2 



Mpc, 



(2) 



where c s is the sound speed, p is the mean density of the 
Universe, f2 m is the matter density parameter, T is the mean 
IGM temperature and h is the Hubble constant of units 
of 100km s~ 1 Mpc~ 1 . However, in the non-linear regime gas 
can collapse to scales smaller than this and the Jeans scale 
becomes a somewhat ambiguous quantity. The effective non- 
linear Jeans length, X 3 , is defined as the width of a kernel 
of the form [1 + {kX 1 /2ir) 2 ]~ 2 , such that the rms fluctuation 
amplitude of p DM (x) is the same as that of the unsmoothed 
dark matter density filtered with this kernel (see section 4 
for details). On scales larger than the effective Jeans length, 
equation can be written as 



t(z, w) = A(z) 



Pdm(^ x ) 
Pdm 



H[w—x—v p (x),b(x)] dx, (3) 



with 
A(z) 



H(z] 



nm 



0.61 



SOOkms^Mpc" 1 ^ f fl b h 



H(z) 



0.02 



/ r p hot \ 

Vio- 12 s-v 



T 



1.5 x 10 4 K 



1+2 



(4) 

where Sl b is the baryonic density in terms of the critical 
density, H(z) is the Hubble parameter, and r p hot is the 
photoionization rate per hydrogen atom. The Doppler pa- 
rameter in the last equation depends on nm as 6 oc n^ c 
NH99 defined the local optical depth as, 



-a/2 



f(x) = A 



p{x) 



(5) 



which is related to the observed optical depth r by a convo- 
lution with a Voigt profile as described in equation Q. 

NH99 have presented a direct Lucy-type iterative 
scheme (Lucy 1970) to recover the optical depth and the 
corresponding mass and velocity fields in the LOS from the 
normalized flux, F = exp(— r). In our tests with hydro- 
simulation (described in more detail in S I4.6L we found that 
estimating the velocity field from the spectra itself is not 
very accurate. We will therefore be less ambitious here and 
use the algorithm of NH99 to recover the gas density in red- 
shift space. We will later use hydrodynamical simulations to 
address the effect of redshift space distortions (see § 14.4ft 
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Table X. A table describing the redshift bins used in this paper. 
Columns number 1,2,3 and 4 show the redshift bin number; its 
mean, minimum and maximum redshifts, respectively. Columns 
5 and 6 show the number of spectra included in each bin and its 
total length in kms . 



z bin 


<z> 


z min 


z max 


# of Spectra 


Length in km/s 


1 


3.29 


3.0 


3.6 


7 


202658 


2 


2.85 


2.6 


3.1 


10 


308166 


3 


2.48 


2.2 


2.7 


10 


308002 


1 


1.96 


1.6 


2.2 


13 


486124 



NH99 showed that the density field can be successfully 
recovered below a threshold value above which the corre- 
sponding flux saturates. The NH99 reconstruction method, 
therefore, imposes an effective upper limit on the recovered 
optical depth in these regions. This will inevitably affect the 
amplitude of the measured power spectrum but as will be 
shown later, not its shape. We use this to determine the 
shape and normalization of the power spectrum separately 
in two steps. 



3.2 The shape of the ID power spectrum of the 
gas density 

Due to the 1-dimensional nature of the Lya forest, the mea- 
sured power spectrum at a given wave-number, fc, is neces- 
sarily the line of sight power spectrum, P 1D (fc) . The relation 
between the ID power spectrum and the three dimensional 
power spectrum, P 3D (k) is (Kaiser & Peacock 1991), 



P 3D (fc) 



2tt dP 1 



k dfc • < 6 > 

The useful dynamical range covered by the Lya forest 
data is restricted to 0.1Mpc _1 h < k < 10Mpc _1 h. The lower 
limit comes from the limited length of the observed QSO 
spectra. The upper limit is imposed by a combination of 
the effective Jeans scales below which the pressure gradients 
wipe out the small scales fluctuations in the baryons and the 
contamination induced by metal lines. In order to measure 
the ID power spectrum from the 31 QSO spectra, we have 
divided the Lya data into four redshift bins (see tableland 
figure 0. For each redshift bin, the sections that are taken 
into account are those that belong to Lya spectra that have 
more than 40% of their total length in the redshift bin at 
hand (see table 0. In order to measure the power spectrum 
within a certain redshift bin, we first calculate the local op- 
tical depth for each section within this bin using the NH99 
method. For a given value of a, the ID matter density is 
calculated from equation The recovered density section is 
then Fourier-transformed and the power spectrum at a given 
wavenumber is obtained. We have estimated the mean and 
the measurement error of the ID power spectrum at a given 
wavenumber from all the individual spectra within the bin. 
The variance is calculated using two independent methods 
that give very similar results. The first is a standard devia- 
tion measurement of all the power spectrum values at a given 
wavenumber while the second uses a bootstrap technique. 

The recovered optical depth f is a good approxima- 
tion to the true field only in regions with f smaller than a 
certain value. For large optical depths, the recovered f typ- 



ically underestimates the true field. As in NH00, we define 
a truncated local optical depth ft as ft = f for f < f c , 
and ft = otherwise. In order to test how this cut-off in 
the optical depth affects the shape of the ID power spec- 
trum, the power spectrum inferred from a truncated optical 
depth distribution is measured for a range of cut off values. 
The resulting dimensionless ID power spectra 1 are shown 
in Figure|5| all renormalized to the same fiducial amplitude. 
The left panel is for (z) = 2.48 and the right panel is for 
(z) — 2.85. The error bars show the variance around the 
mean at each point. Notice, that these errors are indepen- 
dent. Changing the cut-off value of the optical depth has 
very little effect at fc 10Mpc _1 h . This is perhaps not too 
surprising as the regions in the spectrum where the flux is 
saturated have a small volume filling factor. 

At small scales with fc ^ lOMpc - h metal absorption 
lines within the Lya -forest contaminate the signal (e.g. Kim 
et al. 2004). We therefore truncate the recovered ID power 
spectrum at this wavenumber. 



3.3 The amplitude of the ID power spectrum of 
the gas density 

Most of the regions where the Lya flux is saturated corre- 
spond to scales comparable to the Jeans scale. These re- 
gions are scattered roughly randomly across the spectrum. 
The amplitude of the ID power spectrum of the gas density 
can thus not be reliably measured directly from the optical 
depth truncated at a certain cut-off value. 

We will use instead the probability density function 
(PDF) of the recovered optical depth which is related to 
the PDF of the gas density. We thereby use the fact that 
the shape of the PDF changes with increasing amplitude of 
the power spectrum. NH00 showed that the first two mo- 
ments of the PDF of the gas density can be recovered rather 
well from the moments of the probability distribution of the 
truncated optical depth. The main points are summarized 
in the following. 

NH00 define the moments of the truncated optical 
depth f t which can be written in terms of V, the density 
PDF, as, 



< T t >= A 



{l + S) na V(S)dS, 



(7) 



where 8 = p/p — 1 is the density contrast and S c — 
{f c /A) 1/a -l. Further, v is defined as, v — [\n(p/p)~- p{\/ p,2, 
where p,i and p,2 are the average and rms values of In (1 + S). 
NH99 had shown that the PDF of the DM density smoothed 
on the scale relevant for the Lya forest (the effective Jeans 
scale) can be reasonably well approximated by a log-normal 
distribution (e.g. Bi & Davidsen 1997 or Sheth 1998 and 
Gaztanaga & Croft 1999 for different forms of the PDF). For 
a log-normal density distribution, V{y) = exp(—u 2 /2)/\/2n, 
the truncated moments in (0 can be written as, 



< Tt > 



A n (I 



2 2 2 

n a p,2 + napi 



1 + erf 



nap2 



V2 



1 We will generally use the dimensionless ID and 3D power spec- 
tra to express our results. These are defined as tt~ 

-i fcp iD for the 
ID case, and A 2 = (27r 2 )- 1 fc 3 P 313 for the 3D case. 
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Figure 2. The measured shape of the ID power spectrum of the gas density at two different redshifts for a range of f -cutoff values. The 
four cutoff values are f cuto ff = 1, 2, 3 and 10 and the spectra are normalized to the same amplitude. The power spectra were calculated 
from all sections of the absorption spectra in the specified redshift bin. The error bars show the 1 — a errors. For k ,$ 10Mpc _1 h where 
the Lya forest is dominant, the changes in the shape of the power spectrum shape are very small. At small scales with k ^ 20Mpc~ 1 h, 
however, associated metal absorption results in a clear change in the slope when the cutoff value is varied. Note that only wavenumbers 
with values 10Mpc _1 h are used to estimate the 3D matter power spectrum. 



where v c is the value of v corresponding to 6 C ■ By expressing 
v in terms of f in J8J the truncated moments can be written 

as 

< Tt™ > = — exp ^in 2 a 2 ^2 — + 7iln.4^ x 



■ erf 



In t c — na fi 2 — In A + ■ 



Q/i 2 



V2 



(9) 



Here the relation fii — —$/2, that follows from the condi- 
tion < 5 >— for the log- normal distribution, is used. The 
moments of the truncated optical depth distribution depend 
on four parameters A, a and f c . The parameter f c is 
chosen such that for f < f c the local optical depth does not 
suffer from the biases introduced in saturated regions. As 
apparent from equation @ there are two basic degeneracies 
leaving two independent parameters, 



B = ]nA-at&/2, 



(10) 



NH00 showed that the moments of ft can then be written 
in terms of these parameters as 



1 

< T t >= - exp 



n 2 C 2 



+ nB 



erf 



In f c — nC 



B 



V2C 



The first two moments, < f t > and < f t 2 >, are sufficient to 
determine the parameters B and C. From these one can then 
infer the rms fluctuation amplitude of the gas density aj 
and the normalization constant of the optical depth A. The 
rms fluctuation amplitude of the gas density aj is related 
to the amplitude of the ID power spectrum of the gas by a 
simple integration, 



2 

0"J 



(2tt) 3 



P 3D d 3 k 



P 1D dk 



(12) 



where equation |S| is used to obtain the equality on the right 
hand side. 

The normalization adds another uncertainty to the 
power spectrum calculation. We estimate this uncertainty to 
be w 40% (see NH00) for a spectrum as long as Q1422+231. 
To account for the larger total redshift path in each redshift 
bin compared to that of Q1422+231 we assume that the am- 
plitude errors follow a Poisson distribution, i.e., that they 
scale inversely with the square root of the combined length 
of the spectra in each redshift bin. This scaling reduces the 
error associated with the normalization to ~ 10%-20% for 
each redshift bin. This error is added to the error due to the 
shape measurement. Note however, that unlike the shape 
measurement error, the error due to the normalization is 
added to all data points equally and is highly correlated. 

We have estimated the influence of the truncation at 
k — 10Mpc _1 h on the calculation of aj from equation 1121 
extrapolating the dependence of the power spectrum on k 
at k ^ 2Mpc~ h. At these wavenumbers the ID power spec- 
vH}um scales roughly as fc -3 (see figure[5J. This yields a con- 
tribution of the order of a few percent compared to the con- 
tribution of the sampled range. Note that, since the slope of 
the power spectrum is expected to further steepen at larger 
wavenumbers due to the effective Jeans scale cutoff, this is a 
conservative upper limit. This will be further demonstrated 
with simulated spectra (see right panel of figure [7J. 

Figure [3] shows the measured ID power spectra of the 
gas distribution (solid curves), assuming the a values in- 
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Figure 3. The solid curves shows the measured ID power spectra of the gas density for four different redshifts assuming the a values 
indicated on each panel. The error bars include the combined contributions from the measurement of shape and amplitude of the power 
spectrum to the errors as described in the text. The central dashed curve in each panel shows the best fitting models obtained from the 
likelihood analysis. The other two dashed curves show the ~ 1 — a range of allowed values of its- The parameters of the best fit models 
are annotated on each panel. 



dicated on each panel in the four redshift bins. The error 
bars shown are those associated with the uncertainties in the 
shape and amplitude measurements. The figure also shows 
as the dashed curves the best fitting models obtained from 



the likelihood analysis, which will be described in the fol- 
lowing two sections. 
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4 MODELLING THE ID GAS POWER 
SPECTRUM 

4.1 General Considerations 

Equation|S|can be used in order to obtain the 3D power spec- 
trum from the measured ID power spectrum. Unfortunately 
however, the measured ID power spectrum is already a noisy 
quantity and differentiation will introduce more uncertain- 
ties. This approach is explored in §|5|where we compare the 
3D power spectrum obtained here to the 3D matter power 
spectrum inferred from the flux power spectrum by Viel et 
al. (2004). In order to quantitatively constrain the parame- 
ters describing the matter power spectrum and the thermal 
state of the gas, we use instead a likelihood analysis. The 
likelihood analysis is described in section lo~T1 

To obtain a realistic model for the ID power spectrum 
of the gas density we start with the linear 3D matter power 
spectrum in real space. From this the non-linear 3D power 
spectrum of the gas distribution is obtained by taking into 
account the effects of non-linearity, redshift distortions and 
gas vs. dark-matter bias. The model for the ID power spec- 
trum of the gas in redshift space is then readily obtained by 
using the integral form of equation |S| 

p 1D (fc) = _ / P 3D (k')k'dk'. (13) 
2n Jk 

4.2 The Linear 3D Matter Power Spectrum 

We restrict our analysis to the generalized family of CDM 
cosmological models, allowing variations in the mass density 
and vacuum energy density parameters. The general form of 
the power spectrum of these models is, 

P 3D (k) = ylT 2 (fi,fi m ,« b ,h;fc)ifc™ s . (14) 

The CDM transfer function, proposed by Sugiyama (1995), 
is adopted, 

T(k) = 
V ' 2.Mq 

[l + 3.89q + (I6.U7) 2 + (5.46g) 3 + (6.71g) 4 ] _1 (T5) 

q = k [Qh exp(-fi b - h 1/2 Q h /Q) (/iMpc" 1 )]" 1 , (16) 

where Q, S7 m and fib are the total, mass and baryonic density 
parameters, respectively; h is the Hubble constant in units 
of 100km s^Mpc" 1 and n s is the power law index of the 
primordial power spectrum. The parameter A is the power 
spectrum normalization factor, which will be expressed here 
in terms of the other cosmological parameters and as - the 
rms density fluctuations within top- hat spheres of 8h _1 Mpc 
radius. 

The parameters are varied such that they span a 
range of plausible CDM models. In all cases, the bary- 
onic density is assumed to be f2b = 0.02h~ 2 , the value 
currently favoured by primordial nucleosynthesis analysis 
(e.g., Walker et al. 1991, Buries & Tytler 1998) and the 
WMAP cosmic microwave background (CMB) data (Bennet 
et al. 2003). The value of h is fixed to 0.7 and the Universe 
is assumed to be flat. 

We have focused in our study on a few cosmological 
parameters, f2 m , as and n B . It is important to note that 
some of the free parameters of the power spectrum models 



are degenerate. For example at the range of wavenumber 
constrained by Lya forest data the amplitude of the power 
spectrum is degenerate with the power spectrum power law 
index, n s . We have therefore fixed the value of the power 
law index to unity in most of the analysis of the real data. 
We will discuss the influence of changing the value of n s on 
the results for a few cases. 



4.3 Nonlinear effects 

As discussed earlier, the fluctuations in the Lya optical 
depth follow those of the dark matter down to the Jeans 
scale, which is of the order of 1 h _1 Mpc. At redshifts of 2- 
4, these scales have already entered the quasi-linear regime. 
It is therefore necessary to account for the non-linear evo- 
lution of the density field. Peacock and Dodds (1996) have 
developed a simple recipe for mapping the linear 3D mat- 
ter power spectrum onto the quasi linear regime. The actual 
mapping used here is the one described in Peacock's (1999) 
book which has slightly modified, and more accurate param- 
eters than the one given in Peacock & Dodds (1996). The 
recipe is accurate to few percent for a very large family of 
CDM power spectra. For more details, see Peacock (1999) 
and references therein (see also Smith et al. 2003). 

It is worth noting that it should not be necessary to 
follow the evolution of the power spectrum to the highly 
non-linear regime. High density regions will not give rise to 
the low column density spectral lines we use in our analysis. 
The influence of saturated lines is further diminished by in- 
troducing a severe cutoff to the optical depth (see figure 

4.4 Redshift distortions 

After taking the non- linear effects into account, the theoret- 
ical 3D real-space power spectrum is transformed into red- 
shift space. The "distortions" caused by this transformation 
- normally called redshift distortions - will significantly bias 
the power spectrum measurement. On linear scales, redshift 
distortions cause an enhancement in the measured power 
spectrum with a constant factor that depends on the value 
of Qm at a given redshift (Kaiser 1987). In the highly non- 
linear regime, redshift distortions tend to dilute the distribu- 
tion along the line of sight and create the so called "fingers 
of God" . Redshift distortion in the highly non-linear regime 
thus lead to a suppression of the power spectrum relative to 
the real-space power spectrum. The effect of peculiar veloc- 
ities in this regime can be modelled in a statistical manner 
with a Gaussian fit to the ID pairwise velocity distribution 
(Davis & Peebles 1983). In the case at hand however, most of 
density fluctuations are in the linear to quasi-linear regimes. 
Unfortunately, there is no good analytical description of red- 
shift distortions in the quasi-linear regime. 

To account for the redshift distortions we have thus 
used the hydrodynamical simulations described in section 
§ 14.61 to determine an empirical relation between the real- 
space and redshift power spectra (see figure ■ Note that 
this correction is dynamical and as such is not dependent on 
the details of the gas physics. 

The upper left panel of figure0|shows the dimensionless 
ID power spectra with (dashed curves) and without (solid 
curves) peculiar velocity distortions. The magnitude of the 
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Figure 4. Redshift distortions of the ID and 3D power spectra and their errors as estimated from the numerical simulations. The upper 
left panel shows the real-space (solid curves) and redshift-space (dashed curves) ID power spectra at two redshifts, 2.75 (thin) and 2.25 
(thick). The lower left panel shows the ratio between the dimensionless real space 3D power spectrum and its redshift space counterpart 
for three redshifts, 2.25, 2.75 and 3 as dashed/dotted curves. The solid curves shows the fit to the ratio used in the likelihood analysis. 
The upper right panel shows the errors introduced in the dimensionless 3D redshift-space power spectrum by using this fit. The solid 
curves shows the dimensionless 3D power spectrum at z = 2.75; the dashed curve shows the absolute error and the dashed curve shows 
the relative error which has a maximum value of about 20%. The errors at the other two redshifts are similar. The lower right panel 
shows the actual 3D real space power spectra (solid curves) vs. the reconstructed power spectra (dashed curves) at two redshifts, 2.75 
(thick) and 2.25 (thin). 



effect may appear rather large, however one should keep in 
mind that the ID power spectrum is an integral quantity 
(Eq. 1131 which tends to enhance any systematic biases in 
the 3D power spectrum. 

To avoid this accumulation of bias in the ID power spec- 
trum, the influence of the velocity distortions is modelled 
using the 3D power spectrum. The lower left panel of fig- 
ure 2] shows the ratio between the real-space and redshift 
space 3D power spectra in the relevant range of wavenum- 
bers as deduced from the numerical simulations. The ratio 
is shown for 3 different redshifts. At the largest scales the 
ratio flattens to a constant value that is comparable to the 
linear effect described by Kaiser (1987). On smaller scales 
the trend reverses and the redshift space power spectrum is 
suppressed relative to the real-space power spectrum. We fit 
the ratio between the two with the following simple formula, 

= 0.535 (l + (fc/fc ) 3 )*, (17) 

where P r (k) and P s (k) are the real-space and redshift space 



power spectra respectively; ko = 2Mpc _1 h is the scale rel- 
evant for the transition from linear to quasi-linear regime. 
The fit is shown with the solid curve. From figure^Jit is clear 
that the required correction for the ratio is almost indepen- 
dent of redshift. A correction with the same functional fit 
can thus be applied to all redshifts. The errors introduced by 
the fitting formula are shown in the upper right panel. The 
relative error amounts to a maximum of about 20%. The 
lower right panel shows how well the correction recovers the 
real 3D redshift space power spectrum. 

For flat cosmological models, the Universe at the red- 
shifts of the Lya spectra is very similar to an Einstein-de- 
Sitter model. The correction proposed in equation 1171 will, 
therefore, not be very sensitive to the exact values of cos- 
mological parameters, possibly with the exception of ag. We 
have shown that the relation holds at three different red- 
shifts. The functional form of equation 1171 should thus be a 
good fit which depends only weakly on cosmological param- 
eters. 
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4.5 Gas vs. Dark-Matter Power Spectra 

The next issue is to model the bias between the gas and dark 
matter distribution. On scales larger than the Jeans scale, 
numerical simulations have shown that pressure effects are, 
as expected, negligible. The gas faithfully traces the dark 
matter on these scales. On scales comparable or smaller than 
the Jeans scale, however, the distribution of gas deviates 
from that of the dark matter. At these scales, the pressure 
becomes important and prevents the gas from contracting, 
smoothing out all the small scale fluctuations. 

An approximate functional form for the effect of Jeans 
smoothing on the power spectrum can be calculated from 
linear theory (e.g. NH99). This is, however, not adequate 
to describe the quasi-linear regime, relevant for this study. 
We have therefore again used the numerical simulations de- 
scribed in § 14.61 to obtain a more accurate fit to the ratio 
between the power spectra of gas and dark matter densities. 
The solid curves in figure |S] show this ratio for four red- 
shift bins as measured directly from the simulations. Notice 
the slight enhancement at k w 20Mpc~ h. The long-dashed 
curves are the best fit functions based on linear theory, which 
do not reproduce the enhancement. The dashed curves show 
the results of fitting the simulation data with a fitting func- 
tion of the form. 



,313 p dm (1 ■ 



B(z)k 2 



1 + 



' kXj 

, 27T 



(18) 



where PgEi and PxJm are the 3D gas and dark-matter power 
spectra respectively, and Xj is the effective Jeans scale, and 
B(z) = 0.1447-0.0186 z. The function B(z) is later used in 
the likelihood analysis where the value of X.j was kept as a 
free parameter. 

It worth noting that we have used simulations with and 
without cooling. The main difference we found is that sim- 
ulations without cooling require the addition of the term, 
B(z) k 2 , to the fit. Whereas the simulations with cooling re- 
quire no such addition to the normal Jeans cutoff function. 
In addition to that, we suspect that the SPH simulations 
with cooling-functions suffer from excess cooling of the gas, 
a well known issue in SPH simulations. For those reasons, 
we adopt the more general fit that the hydrodynamical sim- 
ulations without gas cooling produce. 



4.6 Testing the recovery of the ID power 
spectrum from the optical depth with 
numerical simulations 

In order to test how well our likelihood analysis recovers the 
non-linear ID gas power spectrum in redshift space we have 
tested it with artificial spectra produced from state-of-the 
art SPH simulations. A suite of simulations with varying 
particle numbers, resolution and boxsize have been carried 
out with the parallel TreeSPH code GADGET-2 (Springel, 
Yoshida & White, 2001; Springel 2005). GADGET-2 was 
used in its TreePM mode which speeds up the calculation 
of long-range gravitational forces considerably. The simu- 
lations were performed with periodic boundary conditions 
with an equal number of dark matter and gas particles 
and used the conservative 'entropy-formulation' of SPH pro- 
posed by Springel & Hernquist (2002). The mean UV back- 



ground produced by quasars as calculated by Haardt & 
Madau (1996) has been assumed. This leads to reionization 
of the Universe at z ~ 6. The simulations were run with 
the equilibrium solver for the thermal and ionization state 
implemented in GADGET-2. The heating rates at z > 3.2 
were increased by a factor of 3.3 in order to achieve temper- 
atures which are close to observed temperatures (Schaye et 
al. 2000, Theuns et al. 2002, Ricotti, Gnedin & Shull 2000). 
At 2 < 6 the power-law index of the gas density temperature 
relation is 7 ~ 1.6, where T = T (1 + <5) 7(z)_1 . 

To maximise the speed of the simulations, a simplified 
star-formation criterion in the majority of the runs is em- 
ployed. All gas at densities larger than 1000 times the mean 
density was turned into collisionless stars. The absorption 
systems producing the Lya forest have small overdensity so 
this criterion has little effect on flux statistics, while speed- 
ing up the calculation by a factor of ~ 6. All feedback op- 
tions of GADGET-2 in the simulations have been turned off. 

We have run four simulations with a box size of 60, 30, 
15 and 10 comoving Mpc, respectively. The three larger 
simulations were run with 2 x 400 3 particles including gas 
cooling. The simulation with a box size of 10 comoving h~ 
Mpc was run with 2 x 200 3 particles and without radiative 
cooling. This was done in the smallest box size simulations 
only in order to address the effect of the Jeans smoothing 
without allowing the gas to radiatively cool. In these sim- 
ulations the thermal state of the gas is set by the equilib- 
rium between photo-heating and adiabatic cooling caused 
by the expansion of the Universe. The cosmological parame- 
ters were chosen to be consistent with the values obtained by 
the WMAP team in their analysis of WMAP and other data 
(Spergel et al. 2003), Q m = 0.26, Ov = 0.74, fi b = 0.0463 
and H = 72 km s" 1 Mpc -1 . The CDM transfer functions of 
all models have been taken from Eisenstein & Hu (1999). 

The left panel of Fig. |7| shows the 3D gas power spec- 
trum for the simulation with a box size of 60 (dotted), 30 
(dashed), 15 (dotted-dashed), and 10 (triple-dotted-dashed) 
h~ x Mpc. Note that the hydro-simulation have insufficient 
dynamic range to simulate the forest down to the Jeans scale 
and capture at the same time the largest structure probed 
by the observed Lya forest. Unfortunately, such a dynami- 
cal range is currently inaccessible with a single simulation. 
We have therefore combined the 3D power spectra of the 
four simulations to obtain an estimate of a power spectrum 
that spans scales that are comparable to those probed by the 
data. Note that this approach can not be applied directly to 
the ID power spectra of the simulations due to the inte- 
gral nature of the ID power spectrum. The best estimate of 
the ID power spectrum expected for infinite dynamic range 
is, thus, obtained with equation 1131 using the combined 3D 
power spectrum of the gas obtained from the simulations. 
This mock ID power spectrum is calculated from the simu- 
lation outputs at z — 2.75. We have added errors that are 
comparable to those found in the data. 

The right panel of Fig. |7| shows the dimensionless ID 
power spectrum obtained by applying an artificial Jeans- 
like cutoff with Xj = 0.8h _1 Mpc to the ID power spectrum 
obtained from the combined 3D power spectrum from the 
simulations (diamond symbols with error bars) . The dashed 
curve shows the non-linear analytical model for the same 
cosmological parameters and amplitude of the linear mat- 
ter power spectrum as implemented in the hydro-simulation 
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Figure 5. Gas vs. dark matter power spectra. The bias between the gas and the dark matter due to Jeans smoothing is shown. The 
solid curve shows the ratio of the gas and dark-matter 3D real space power spectra. The short dashed curve is the fitting function used 
in the current analysis. The long-dashed curve is the normal Jeans scale cutoff function. Note that the hydro-simulations used here do 
not include gas cooling (see text). 



with a Jeans cutoff at Xj = 0.8h~ Mpc. The agreement 
between the non-linear ID power spectrum model of the 
gas density and our best estimate from the set of hydro- 
simulations is, perhaps not too surprisingly, excellent. 

4.7 Summary of the analysis steps 

To summarize the analysis steps: First, the non-linear ID 
power spectrum of the gas distribution is modelled analyti- 
cally. In this way we can easily vary the shape and amplitude 
of the underlying linear matter power spectrum. The ther- 
mal state of the gas is modelled by a thermal smoothing 
function which has one free parameter, the effective Jeans 
length. The shape of this function is derived from our hydro- 
dynamical simulations. Second, the model ID power spectra 
are compared with those recovered from the truncated op- 
tical depth of the Lycv forest data using a maximum likeli- 
hood analysis. These steps are summarized in the flow-chart 
shown in Fig. |S] 



5 CONSTRAINING THE PARAMETERS OF 
THE MATTER POWER SPECTRUM AND 
THE THERMAL STATE OF THE GAS 

5.1 The likelihood method 

For this the following likelihood is maximized, 



£ = COnSt. X e ~5l P ob a - P modelJ C l P ob B - P modcl ) y (1Q) 

where Pj^ 8 and P^dci are the observed and model power 
spectrum vectors and C is the error correlation function. 

As mentioned earlier, given the large uncertainties 
in the data and the various degeneracies in the models 
(e.g., between the power law index, n s and erg) we will 
not attempt to simultaneously constrain all free parameters 
on which the power spectrum models depend. Instead, the 
analysis will focus on three or less free parameters at a time 
while keeping the other parameters fixed. 



5.2 Error estimates 

The errors used in the likelihood analysis come from two 
sources, the ID power spectrum measurement and the un- 
certainties in the theoretical model used to fit it. The error 
in the measurement were discussed earlier in § 13.21 These 
were shown to have an uncorrelated contribution due to the 
power spectrum shape determination and a correlated com- 
ponent due to the amplitude determination. 

The main source of error in the model comes from the fit 
used to account for the redshift distortions which introduces 
a error that can be as high as 20%. To simplify the treatment 
we conservatively fix the relative error to be 20%. This error 
is added coherently to all wavenumbers. 

To summarize, the error matrix, C, is defined as follows: 
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Figure 6. A flow-chart summarising the analysis steps as applied to the data and the model. The sections in the paper in which each 
step is explained are indicated at the upper right corner of each "step-box" . 



C — dia(7(A s i lape )+ (^amplitude + e z — distortions) PobsP obs j (20) 

where A s h apc is the error introduced by the power spec- 
trum shape measurement. We performed a covariance anal- 
ysis which showed that, for a given choice of wavenumbers 
k, this error is not correlated across different wavenumbers. 
£ampiitudo is the relative error due to the amplitude measure- 
ment and e z _distortions is the error produced by the correction 
for redshift distortions. Both errors are highly correlated. 

Typical values of A s h ape are shown in figure |2] and 
amount on average to about 50% at each k (see section S 13.21 
for a detailed description). The value of e z -distortions is taken 
to be 20%. The value of e amp iit u de is taken to be 20% for a 
spectrum of the same length as that of Q1422+231 and is 
assumed to scale with the inverse of the square-root of the 
total length of the spectrum used in each redshift bin. This 
scaling of the errors assumes a Poisson distribution of the 
Lya absorption features. 



5.3 Testing the likelihood analysis with the 
hydro-simulations 

In order to test how well our procedure of constraining the 
model parameter works, we have performed a likelihood 
analysis for a sample of artificial spectra obtained from the 



hydrodynamical simulations. The numerical simulation was 
for a fiat ACDM model with fi m = 0.26, <j 8 = 0.85, n s = 1, 
and has a temperature-density relation that corresponds to 
a — 1.56. For the likelihood analysis, the Universe is as- 
sumed to be flat and that a = 1.56. This leaves f2 m , as and 
the effective Jeans length Xj as free parameters. 

For the analysis of the hydro- simulation, the "compos- 
ite" non-linear ID power spectrum, shown in right panel 
of Fig. |7| is used. Figure |H| shows three two dimensional 
likelihood contours (with marginalization over the third pa- 
rameters) and the one dimensional likelihood (normalized 
to have a maximum of unity) for each of the parameters 
with the other two marginalized over. The top left panel 
shows flm vs. ag. The actual values are denoted by a cross. 
For the correct value of £2 m the actual value of as falls well 
into the 1-cr limits of the recovered value suggesting that the 
method works fine. Note, however, that there is a degener- 
acy of the inferred as with the assumed Q m . Not surprisingly 
the Lya forest data alone cannot constrain both parameters. 
The top middle and right panel of Figure |5] show the like- 
lihood contours of the effective Jeans length X.j vs as and 
Sl m , respectively. Due to the cut-off of the power spectrum 
at small scales, Xj is well constrained. 

The solid curve at the right panel of Fig. |7| shows the 
non-linear model power spectrum for the best fit parameters 
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Figure 7. Left panel: shows the 3D power spectra of the gas density as measured from the simulations with a box size of 60 (dotted), 
30 (dashed), 15 (dotted-dashed), and 10 (triple-dotted-dashed) h~ l Mpc boxsizes. Clearly, none of the simulations spans the dynamical 
range needed to capture simultaneously the large and small scales needed to model the ID power spectrum. The grey solid curve is the 
composite 3D power spectrum which is our best guess of what the non-linear 3D power should look like with infinite dynamic range. 
Right panel: the diamond symbols show the ID power spectrum of the gas density corresponding to the composite 3D power spectrum 
in the left panel (with error bars) with a Jeans cutoff at Xj = 0.8h _1 Mpc. The solid curves shows the best fit power spectrum found by 
the likelihood analysis with cosmological parameters as annotated on the figure. The dashed curve shows the ID power spectrum with 
the actual parameters of the hydro-simulations with Xj = 0.8h -1 Mpc. 



obtained with the likelihood analysis (see The agree- 
ment with the best estimate from the hydro-simulations is 
again very good. The recipes used to model the non-linear 
evolution of the matter power spectrum, the bias between 
dark matter and gas density and the recovery of the pa- 
rameters describing the optical depth distribution, appear 
to work very well. 

We have also tested how sensitive our analysis is to the 
functional form which we use to model the Jeans smoothing 
by setting B(z) — 0. We found no significant difference with 
regard to the deduced cosmological parameters except on the 
value of Xj itself (for more discussion on this see section 0. 
The reason is that even with the canonical functional form 
for the Jean's smoothing deviations from the actual ratio are 
« 15% within the observed k range. This is small compared 
to other uncertainties in the modelling. 

5.4 Constraints on erg, Xj and fi m for a flat 
cosmology 

We now turn to a likelihood analysis of the real data. One 
difficulty is that we will have to assume a value of a, the 
power law index used in equation [S] which relates the local 
optical depth with the underlying density. However as we 
will see later, by comparing results of different redshifts, 
we can actually constrain the evolution of a and thus the 
thermal evolution of the IGM. As discussed earlier, a ranges 
between 1.56 — 2, with the lower limit corresponding to the 
balance between adiabatic cooling and photoheating and the 
upper limit to an isothermal IGM. During a phase in which 



the IGM is rapidly heated (e.g., during the reionization the 
value of a becomes closer to 2). 

Equation |S] suggests that for a given value of the local 
optical depth, f , the amplitude of the density will be larger 
for smaller a and vice versa. In this subsection, a is chosen 
to be either 1.56 or 1.8, for all redshift bins. Also shown are 
results for a 'mixed' case in which a evolves with redshift as 
expected for a scenario in which Hen is ionised to Hem at 
z « 3.2. 

As discussed above, the value of n s is degenerate with 
erg, and we have assumed n a to be fixed to a value of unity. 
We have furthermore restricted our analysis to models with 
a flat cosmology (i.e., Q m + £7a = 1). The free parameters 
of this analysis are S7 m , Xj and erg. 

Figures [5] and 1101 show the results of our analysis as- 
suming q = 1.56 at z = 2.48 and z — 3.29 respectively. The 
same three two-dimensional likelihood contours and the one 
dimensional likelihood curves as in the test with numerical 
simulations are shown in fig. |S] The contours look gratify- 
ingly similar. As in the case of the hydro-simulations there is 
a degeneracy of the inferred erg with the assumed f2 m while 
Xj is tightly constrained. 

The degeneracy between Q m and erg is somewhat weaker 
and the constraints on erg (marginalized over the other two 
parameter) are somewhat stronger than in the analysis of 
the mock spectra. Note that the constraints on Xj are also 
tighter. These differences can be traced back to the shape 
of the thermal cut-off of the non-linear ID power spectrum 
of the gas distribution which is still affected by the limited 
dynamical range of the simulations. 
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Figure 8. The likelihood results for the power spectrum from the hydrodynamical simulations with Q m , o"8 and Xj as free parameters. 
The upper three figure show the 1, 2 and 3-cr likelihood contours for each pair of the parameters after marginalizing over the third. The 
lower three panels show the likelihood function for each of the three parameters (after marginalizing over the other two). The likelihood 
is normalized to have a maximum value of unity. The value of fi m is poorly constrained, but the values of erg and and Xj are well 
constrained and agree with the simulations parameters. 



Figure [TTI summarizes the results for the marginalized 
likelihood of as for all four redshift bins. The upper left 
panel shows the likelihood plots for the 4 bins, calculated 
with a = 1.56. The results for the bins with z =1.96, 2.48, 
2.85, and 3.29 are shown as dotted, dashed, dotted-dashed 
and long-dashed lines respectively; the solid curve is the joint 
likelihood curve. The likelihood lines shift slightly between 
different redshift bins. The upper right panel is the same as 
the upper left panel but for a = 1.8. The preferred value 
of erg drops here, as expected, by about 30%. The lower left 
panel is the same as the previous two panels except that 
a =1.56, 1.7, 1.8 and 1.65 for z =1.96, 2.48, 2.85, and 3.29, 
respectively. Such an evolution is consistent with reioniza- 
tion of Hell at z ~ 3.2. 

The lower right panel of figure fTTI shows the likelihood 
curves of the effective Jeans length, Xj. Here, we assumed 
the mixed a model, i.e., a values as in the lower left panel, 
for each of the four redshift bins. There is a clear evolution 
between the redshift bins with the exception of the z = 2.85 
and z = 2.48 redshift bins where the evolution is very small. 
The effective Jeans length increases with decreasing redshift. 



This is consistent with an increased heating rate at redshift 
> 3.2 due to helium not yet being fully ionized. The in- 
creased energy injection during Hell reionization should lead 
to an expansion of the sheets and filaments responsible for 
the Lya forest. Note that the timescale for Xj to change 
should be a fair fraction of the Hubble time and that Xj 
is therefore not necessarily a good measure of the instan- 
taneous temperature. We have also examined the effect of 
changing n s . Lowering n s to 0.95, results in an increase of 
the estimated power spectrum amplitude as by 4 %. 



In order to show that the models with the preferred 
parameters of our likelihood analysis (discussed later in this 
section) fit the data well, the measured ID power spectra 
shown in Figure[3]are compared with the theoretical models. 
The models are shown as the dashed curves and have three 
erg values, 0.85, 0.95 and 1.05. Notice that f2 m = 0.3. The 
values of a are consistent with the evolution expected for a 
reionization of Hell at z ~ 3.2 (see discussion later in the 
section) . 
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a =1.56 <z> = 2.48 z min = 2.20 z max = 2.70 




0.2 0.4 0.6 0.8 0.6 0.8 1 .0 12 3 4 



Figure 9. The likelihood results for the observed spectra for the redshift bin at z = 2.48 with Q m , erg and Xj as free parameters. 
The assumed value of a is 1.56. The upper three figures show the 1, 2 and 3-<r likelihood contours for each pair of the parameters after 
marginalizing over the third. The value of n m is poorly constrained, but the values of erg and and Xj are well constrained. 



5.5 Constraints on a and Xj with a fixed 
cosmology and normalization 

Does the data prefer certain values of oil In order to an- 
swer this question more definitely, we carried out an analy- 
sis where a and Xj are left free while the cosmological pa- 
rameters n m , £Ia are fixed to 0.3 and 0.7 respectively. The 
amplitude of power spectrum was set to match that of the 
2dF galaxy survey normalization (Cole et al. 2005). Since 
the mean separation between galaxies in the survey is about 
10h~ Mpc, as is not directly measured by the 2dF data. 
To avoid extrapolating to smaller scales, the 030 amplitude 
is used. This amplitude is defined as the rms fluctuations 
within 30h _1 Mpc spheres, a scale directly probed by 2dF. 
The value of 030 measured by 2dF is 0.233 (Percival et al 
2002). The mass-galaxy bias ratio is assumed to be unity 
and any possible errors of the 2DF measurement (which are 
relatively small) is neglected in the likelihood analysis. 

Figure [T2"l shows the likelihood contours for a and Xj 
for each redshift bin together with the marginalized like- 
lihoods. The data prefers values of a within the assumed 
physical limits of 1.56 and 2, and appears to suggest that 
q evolves with redshift as expected if Hen reionization in- 



deed occurred at z ~ 3.2. The inferred values of Xj are 
consistent with the results obtained in the previous subsec- 
tion, and are only weakly dependent on the values of the 
cosmological parameters. 

5.6 Joint constraints of f2 m — erg from Lya and 
WMAP data 

As apparent from Figs. 8 and 9 and discussed in section 
15.41 the fluctuation amplitude of the linear matter power 
spectrum inferred from the likelihood analysis is somewhat 
degenerate with the inferred value for Q m . This degeneracy 
is shown in the left panel of Fig. ll4l which shows the 2-<r con- 
tours for four redshift bins separately. The likelihood con- 
tours of the four bins are in good agreement. We have here 
assumed again the 'mixed' a model for the optical depth 
density relation and made the same assumptions as in sec- 
tion geu 

As shown in the middle panel of Fig. 1141 this degen- 
eracy is orthogonal to a similar degeneracy for the CMB 
data. The dashed curves show the constraints from the 
WMAP data combined with a prior on the Hubble constant, 
H = 72 ± Skms^Mpc" 1 (Freedman et al. 2001). The er- 
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a =1.56 <z> = 3.29 z min = 3.00 z max = 3.60 




0.2 0.4 0.6 0.8 0.6 0.8 1 .0 12 3 4 



Figure 10. The same as figure |3 but for the redshift bin at z = 3.29. The preferred value of erg and Xj here are lower than those of 
figure [5] suggesting a change of a. 



ror contours were calculated using COSMOMC (Lewis & 
Bridle 2002). The curvature of the Universe, fife, was as- 
sumed to lie between —0.3 and +0.3. Without the prior on 
the Hubble constant the contours broaden somewhat but 
show a similar degeneracy. The right panel shows the much 
tighter constraints on as and S7 m obtained by combining 
our Lya forest data with the WMAP and HST key project 
data. The joint analysis yields the values as = 0.9±0.06 and 
fl m = 0.3 ± 0.05 (in good agreement with Viel, Haehnelt & 
Springel 2004, Seljak et al. 2005 and Viel & Haehnelt 2005). 



6 THE THREE DIMENSIONAL POWER 
SPECTRUM 

In principle it would be more convenient to infer the 3D 
matter power spectrum directly from the data. This would 
also facilitate a more direct comparison with the results of 
Viel et al. (2004) who inferred the 3D matter power spec- 
trum from the flux power spectrum using an effective bias 
method calibrated with numerical simulations. However, as 
discussed above, this requires taking the derivative of noisy 
data. We deal here with this problem by assuming that the 
ID power spectrum is an analytic function. We use a generic 



curve to fit the ID data points. The derivative is then easily 
obtained. We have fitted the following functional form to 
the measured ID power-spectrum, 



mi*) 



(21) 



l + {k/ko)f> ' 

where = (Ao, 7, fco, 0) is the free parameters vector. 

The main issue that remains is how to assign errors to 
the derived 3D power spectrum. In order to estimate the 
errors, let 8f(k; fj,) be the uncertainty in the functional fit 
due the errors in the data. Assuming that the fit has the 
correct functional form, this error is due to the uncertainty 
in determining the free parameters, namely, 



(22) 



The value of can easily be estimated from the fitting 
procedure used to obtain the most likely values of no. Here 
we chose to do this with a minimum \ 2 analysis. Once S)j, is 
known it is straightforward to show that the uncertainty of 
the derivative is given by: 



(§£)-*> 



d (df/dm 

dF~ 



(23) 
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Figure 11. The two upper panels show the likelihood distribution for erg with marginalization over £7 m and Xj for different redshifts 
with a = 1.56 (left panel) and a = 1.8 (right panel). The lower panel show the likelihood distribution for erg and Xj for a model in 
which a evolves with redshift from a = 1.65 at (z) = 3.29 rising to a = 1.8 at (z) = 2.85 and then a decreasing to 1.65 and 1.56 at 
(z) = 2.48 and 1.96 respectively. Such an evolution would be expected if He became fully reionized at z ~ 3.2. The key for the line styles 
is given in the upper left panel. 



The right hand side of equation 1231 is readily calculated and 
gives the errors associated with the inverted quantity. 

The ID power spectra are obtained assuming the mixed 
a model for the optical depth density relation. Figure 1131 
shows the functional fit to the observed ID power spectra at 
each redshift bin for the mixed a model. Equations l21l and l23l 
are used to obtain the 3D power spectra. These have then 
be scaled to redshift zero using the linear growth factor and 
corrections for the effects of redshift distortions and Jeans 
smoothing have been made. For the latter the most likely 
Jeans scale found in the likelihood analysis has been used. 



The left panel of figure 1151 shows the 3D power spec- 
tra for the four redshift bins assuming that a evolves with 
redshift as suggested by our likelihood analysis. The figure 
clearly shows that the 3D power spectra from the four red- 
shift bins is consistent down to the scales where the inversion 
becomes unstable and power spectra start to diverge. This 
instability is caused by the insufficient information content 
in the data below the Jeans smoothing scale. As expected 
the scale of instability varies with redshift. Note that the 
wavenumbers k shown here is the linear k, the actual mea- 
sured k is larger by roughly a factor of 5 (see Peacock & 
Dodds 1991). 
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Figure 12. Constraints on a and Xj with fixed cosmological parameters, f2 m = 0.3, Q\ = 0.7 and amplitude of the matter 
spectrum 0-30 as measured by 2dFGRS. 



power 



The right panel shows the average 3D power spectrum 
with errors that reflect the uncertainty in a and the differ- 
ences between the four redshift bins. The uncertainty in a 
is added by taking the amplitude variation within the range 
1.56 < a ^ 1.8. In the averaging procedure, the points at 
scales smaller than the scale of convergence are ignored. The 



right panel also shows the 3D power spectra inferred by Viel 
ct al. (2004; dashed line) from the LUQAS sample and from 
the Croft et al. (2002) sample (dotted line) as reanalysed by 
Viel et al. (2004). There is excellent agreement between the 
inferred linear 3D matter power spectrum of this study and 
those obtained by Viel et al. (2004). 
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Figure 13. The measured ID power spectrum of the gas density in four redshift bins. The dashed curves show the best-fitting models. 
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Figure 14. Left panel: The likelihood contours in the S7 m — as plane, after marginalizing over Xj for four different redshifts for the 'mixed' 
a model discussed in the text. Middle panel: The joint likelihood of the four redshift bins together, with solid contours representing the 
1,2, and 3-cr likelihood. Also shown are the 1 and 2-cr contours (dashed curves) as obtained from the WMAP data, assuming a Universe 
with fij; between -0.3 and 0.3, and a Hubble constant of 72 ± 8kms Mpc . Right panel: The joint likelihood of the Lya data and 
WMAP data (again assuming the Universe to be flat) and an HST prior on the Hubble constant. 
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Figure 15. Left panel: The linear 3D matter power spectrum at z = inferred from observed data at four different redshifts. The power 
spectra have been mapped to z = assuming a flat Universe with f2 m = 0.3 and = 0.7. The power spectra have been corrected for 
redshift distortions, Jeans smoothing and for non-linear effects as described in the text. Notice that the inversion becomes unstable at 
scales smaller than the Jeans scale relevant to the redshift at hand. The right panel shows the mean and error after combining the results 
from the four redshifts (solid curve). For comparison we show the linear matter power spectrum power spectrum inferred by Viel et al. 
2004 from the LUQAS sample (dashed curve) and from the Croft et al. 2002 sample (dotted curve), as re-analysed by Viel et al. (2004). 



Table 2. A table showing the 3D power spectrum as estimated 
from the direct inversion method. The first column shows the 
wavenumber k. The second column shows the 3D power spectrum 
and the third and fourth columns show the lower and upper error 
bars around the mean. 



k 

[Mpc- 1 ^ 


pSD 

[h- 3 Mpc 3 ] 


_ AP 3D 
[h" 3 Mpc 3 ] 


+AP 3D 
[h~ 3 Mpc 3 ] 


0.200000 


3023.49 


1005.55 


1716.95 


0.248189 


1325.90 


577.095 


760.660 


0.307990 


877.984 


372.994 


495.885 


0.382199 


586.556 


210.709 


296.521 


0.474288 


410.985 


129.560 


187.204 


0.588566 


263.541 


84.0443 


114.476 


0.730379 


163.626 


54.6102 


68.7632 


0.906361 


116.986 


41.3544 


48.0453 


1.12475 


70.5957 


24.1203 


27.1044 


1.39575 


52.7193 


16.9332 


18.2188 


1.73205 


35.7962 


9.52174 


10.8408 


2.14938 


24.4715 


3.33336 


5.43795 


2.66727 


15.7595 


2.04052 


3.73897 


3.30994 


10.1860 


1.78284 


3.32257 


4.10746 


6.36426 


5.19802 


8.69673 



This is gratifying as our inversion method does not re- 
quire the assumption of an effective optical depth (see Lidz 
et al. 2005 who come to similar conclusions by combining 1 
and 2pt statistics of the flux distribution). Table [(^tabulates 
our inferred 3D power spectrum and the 1 — a errors. 



7 SUMMARY & CONCLUSIONS 

We have presented a new method for measuring the mat- 
ter power spectrum from Lya forest data. The ID density 
field is obtained from the reconstructed line-of-sight optical 
depth for Lya absorption, assuming a simple power law re- 
lation between density and optical depth. We thereby follow 
NH99 & NH00 and introduce a cut-off in the optical depth 
to handle the saturation effects of the flux in high density 
regions. The shape of the non-linear ID power spectrum of 
the gas density does not depend on the value of the optical 
depth cut-off. This allows us to the derive the amplitude and 
the shape of the ID power spectrum of the gas density sep- 
arately from 2pt and lpt statistics of the flux. The shape of 
the power spectrum is calculated directly form the recovered 
line-of-sight gas density while the amplitude of the power 
spectrum is derived from the width of the PDF of the gas 
distribution which we assume to follow a log-normal distri- 
bution. In this way we can measure shape and amplitude of 
the ID power spectrum of the gas density without the need 
for determining or assuming a mean flux level, which has 
proven problematic in previous measurements of the mat- 
ter power spectrum utilizing the flux power spectrum. Note 
that the inferred amplitude of the power spectrum still de- 
pends on the assumed power law index a, which describes 
the relation between the neutral and the total gas density 
and depends on the thermal state of the gas. 

We have then compared the inferred non-linear ID 
power spectrum of the gas density to an analytical model 
of the gas power spectrum using a likelihood analysis. We 
have thereby used state-of-the art hydrodynamical simula- 
tions to model the redshift-distortions and the bias between 
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gas and dark matter distribution and to test our method. 
The Lya -forest data cannot constrain all relevant parame- 
ters simultaneously. We have therefore restricted our likeli- 
hood analysis to flat cosmological models. We only vary a 
subset of parameters at any given time setting the remain- 
der to plausible values. We have also performed constraints 
from a joint analysis of our results from the Lya forest data 
and other data (2dF,CMB,HST). Finally we have obtained 
a direct estimate of the 3D matter power spectrum from the 
ID gas power spectrum and compare it to previous estimates 
from Lya forest data utilizing the flux power spectrum. 
Our main results can be summarized as follows. 

• By enforcing the cosmological parameters to be the 
same in all four redshift bins of our likelihood analysis, we 
found evidence for evolution of the density-temperature re- 
lation and the thermal smoothing length of the gas. The 
inferred evolution is consistent with that expected if the 
reionization of Hen occurred at z ~ 3.2. 

• Assuming that the Universe is flat and assuming a fixed 
value for n s = 1, we find that the fluctuation amplitude of 
the matter power spectrum, erg = 0.9(f2 m /0.3) -0 ' 3 where the 
value of a changes as a function of redshift. The la error on 
as at fixed fi m and a is about 10%. The thermal smoothing 
length Xj is also found to be tightly constrained (to within 
10 percent). 

• A joint analysis of the Lya forest and the WMAP CMB 
data together with a prior on the Hubble constant, yields 
tight constraints on the fluctuation amplitude and the mat- 
ter density, a 8 = 0.9 ± 0.09, O m = 0.30 ± 0.05. 

• The inferred linear 3D matter power spectrum agrees 
well with that obtained by Viel et al. 2004 with a very dif- 
ferent analysis technique. 

The good agreement of our results for the amplitude 
and shape of the matter power spectrum with those of pre- 
vious studies of Lya forest data is very reassuring, as the 
systematic uncertainties differ significantly for the different 
methods employed by these studies. The independent con- 
straints on the thermal state of the gas suggest that the 
inferred peak in the photo-heating rate of helium at z ~ 3 
has affected the flux distribution of the Lya forest in a mea- 
surable way. 
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